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Abstract: We introduce a notion of viscosity solutions for a nonlinear degenerate diffusion 
equation with a drift potential. We show that our notion of solutions coincide with the weak 
solutions defined via integration by parts. As an application of the viscosity solutions theory, 
we show that the free boundary uniformly converges to the equilibrium as t grows. In the case 
of a convex potential, an exponential rate of free boundary convergence is obtained. 

1 Introduction 

Consider a C 2 function : R n — ► R, and consider a nonnegative, continuous function 

po(x) : R n — ► R which has compact support Qq. In this paper we study the porous medium 
equation with a drift 

Pt = A(p m ) + V-(pV$), (1.1) 
for m > 1, with initial data po(x). Note that, at least formally (and proven in [3] for the 



weak solutions), the solution of (1.1) preserves its L 1 norm. 



It will be convenient to change from the density variable to the pressure variable 

u=—o m -\ u = ^-p™-\ (1.2) 
m—1 m — 1 

so that the equation becomes 

(PME-D) Ut = (m - l)uAu + \Vu\ 2 + Vu • V$ + (m - l)itA$ 

(for more on the density to pressure transform see e.g., the discussions in [5]). We consider 
continuous and nonnegative solutions in the space-time domain Q = R n x (0, T) for some 
T > 0, with prescribed initial conditions u(x, 0) = uo(x) G C(R n ). 

When $ = 0, (PME-D) is the widely-studied Porous Medium Equation (PME): We refer 
to the book jT5] for the references. Moreover when V = \x\ 2 , (PME-D) is obtained as a 
re-scaled version of the (PME) with new variables 

0(77, t) := t~ a u(x,t),r] = xt' (3 ,T = Int. 



where u solves (PME). This suggests that the local behavior of (PME-D) is similar to that 
of (PME), with perturbations due to the inhomogeneity of <£. We will illustrate this fact in 
the construction of various barriers in section 2 and 3. 

The weak solution theory for (PME-D) in the case of bounded domains has been developed 
in jl] and [8] . Also, in [11] , existence and uniqueness of solutions are established for the full 
space case under reasonable assumptions (either the initial data is compactly supported or 
the potential has less than quadratic growth at infinity). 

Further, uniform convergence to equilibrium for (PME-D) have also been shown in [3] (see 
Theorem 3.1). In |12| . the connection between the (PME) and the nonlinear Fokker-Planck 
equation is established, which facilitates the use of the entropy method to derive an explicit 
L 1 rate of convergence. In |13j . an extensive study is made of a general form of the nonlinear 
Fokker-Planck equation, i.e., pt = V • (V(^(p) + pV^) with suitable assumptions on <p and 
exponential L 1 rate of convergence is obtained. (PME-D) falls under the framework of |13j . 
and in fact it is the case that almost all of our results would also go through for a general 
equation of this form, but for ease of exposition we will restrict attention to (PME-D). 

We introduce a notion of viscosity solution for the free boundary problem associated with 
this equation, which we will show to be equivalent to the usual notion of weak solutions - 
see [7J for the general theory of viscosity solutions. Note that, formally, the free boundary 
d{u > 0} moves with the outward normal velocity 

V= T ^ 1 = (D U + D<S>).^ = \D U \ + D$ DU 



\Du\ \Du\ \Du\ 

where the first equality is due to the fact that u = on T(n). In this regard we closely 
follow the framework and arguments set out in [6] (see also p3] and [5] ) , where the viscosity 
concept is introduced and studied for the Porous Medium Equation. We point out especially 
that [5] extends the result of [6] to the case where the diffusion term is multiplied by more 
general nonlinearities; our focus, however, is on the added drift term, which introduces 
spatial inhomogeneities. The key utility of the viscosity concept here is that we will be able 
to describe the pointwise behavior of the free boundary evolution by maximum principle 
arguments with local barriers. As an application, we are able to extend the results of jl] 
and [T3] to a stronger notion of free boundary convergence. 

Main Theorem. There exists a viscosity solution of u (PME-D) with f ( m ^uo) 1 ^ m ~ 1 dx = 
f po(x)dx = tuq. Further: 

(a) u is unique and coincide with weak solutions studied in J^J and }13f . 

(b) If \D&\ > except at x = xq where $ achieves its minimum, then from 01/ there 
exists a unique Cq > depending only on thq such that u uniformly converges to 
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(Co — $(x))+. Then we have 



sup d(y, r(uoo)) — > as i — > oo. 

y£d{u(;t)>0} 



(c) For a strictly convex $ with ko < A<£, there exists K,a > depending on uq, ko, the 
C 2 -norm of and n such that 



Remark. 1. We point out that the free boundary convergence may not hold if \D$>\ vanishes 
at some points, even though the uniform convergence of the solution still holds. 

2. In the case of <J>(x) = \x\ 2 (that is for the renormalized (PME)) Lee and Vazquez I16ij 
showed that the interface becomes convex in finite time. It is unknown whether such results 
hold for general convex potentials: we shall investigate this in an upcoming work. 

2 Viscosity Solution 

In this section we introduce the appropriate notion of viscosity solution for (PME-D) and 
show that it is equivalent to the usual notion of weak solution. Our definition descends from 
those in [6] and [H]. For more details we also refer the reader to the definitions, discussions 
and results in j5]. 

2.1 Definition and Basic Properties 

Let Q := M n x (0, oo). For a nonnegative function u(x, t) in Q, we denote the positive phase 



As in [6], to describe the free boundary behavior using comparison arguments we need an 
appropriate class of test functions to deal with the degeneracy of (PME-D). 

Let E be a smooth, cylinder-like domain in W 1 x [0,oo), i.e., 



sup d(y,T(u 00 )) < Ke 



-at 



n(u) = {u> 0}, Q t (u) := {x : u(x, t) > 0} 



and the free boundary 



T(u) = SO(u), T t (u) := dQ t (u). 




E(£) x {t}, where E(£) is a smooth domain in ~EL n . 



(2.1) 



ti<t<t 2 
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Definition 2.1. A nonnegative function u G C 2 ' l ({u > 0} PI S) is a classical free boundary 
subsolution in S if u satisfies (PME-D) with < replacing = in the classical sense in {u > 
0} n S, > on T(n) n S with the outward normal velocity 

F < (3\Vu\ + V$ • on r(«) 

or equivalently, 

< /3|Vw| 2 + V$ • Vu on 

We define a classical free boundary supersolution by replacing < with >. Finally, u is a 
classical free boundary solution if it is both a sub and supersolution. 

Before proceeding further it is convenient to introduce some auxiliary definitions. 

Definition 2.2. Let ip be a continuous, nonnegative function. Now if ip is another such 
function, then we say that ip touches ip from above at (xo,to) in £ if <p — ip has a local 
minimum zero at (xo,to) in S n {t < to}. We have a similar definition for p touching ip 
from below. 

Definition 2.3 (Strictly Separated Data). For two nonnegative functions u,v : W l — > R, 
we write uo -< fo if the following holds: supp(uo) is compact and supp(uo) C Int(supp(«o)) 
and inside supp(iio), uq(x) < vq(x). 

We note that e.g., due to the maximum principle, a classical free boundary subsolution 
that lies below a classical free-boundary supersolution at time t\ > cannot cross the 
supersolution from below at a later time t<i > t±. This observation leads to a notion of 
viscosity solution which takes into account the free boundary: 

Definition 2.4. Let u be a continuous, nonnegative function in Q. 

o u is a viscosity subsolution of (PME-D) if, for any given smooth domain S given in 
(2.1), for every ip G C 2,1 (E) that touches u from above at the point (xo,to), we have 

ipt < apAp + (3\X7p\ 2 + V ■ {pV$) (2.2) 

o ii is a viscosity supersolution of (PME-D) if, for any given smooth domain S given in 
(2.1), 

(i) for every ip G C 2,1 (X!) that touches u from below at the point (a?o,io) G n S, 
we have 

VJ t >a99A(^ + /3|Vv9| 2 + V-(99V<I>) (2.3) 

(ii) for every classical free-boundary subsolution <p in S, the following is true: If 
9? -< n on the parabolic boundary of S, then 99 < u in S. That is, every classical 
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free-boundary subsolution that lies below u at a time t\ > cannot cross u at a later 
time *2 > ii- 



o n is a viscosity solution of (PME-D) with initial data uo if u is both a super- and 
subsolution and u uniformly converges to uq as t — > 0. 

Remark 2.5. /n general one can define viscosity sub- and supersolutions respectively as 
upper- and lower semicontinuous functions. Such a definition turns out to be useful when 
one cannot verify continuity of solutions obtained via various limits. This problem does 
not arise in our investigation here thanks to 141, and our definition assumes continuity of 
solutions. 



It is fairly straightforward to verify that a classical free boundary sub- (super)solution is 
also a viscosity sub- (super)solution. 

Lemma 2.6. If w is a classical free boundary sub- (super) solution to (PME-D), then w 
is also a viscosity sub- (super) solution. 



Proof. We will be brief: The subsolution case presents no difficulty since if contact with 
some cp £ C 2,1 (S) occurs in f2(u>) then we use the fact that w is classical there, whereas 
no contact can occur on the free boundary unless |Vto| = 0, in which case the differential 
inequality is satisfied since then cp = |Vy| = and cp t < 0. 



If w is a classical free boundary supersolution, then (i) in Definition 2.4 follows as before. 
To see (ii), let us note that if cp is a classical free boundary subsolution which crosses 
w, then since the free boundary is C 2 , Hopf's Lemma implies that at the touching point 
\Vcp\ < \Vw\ (see e.g., [H]). O n the other hand, since (p started below w, at the touching 
point we must have v n (cp) > v n (w), which leads to a contradiction since it is also the case 
that we have ^« = V^. 

□ 



Next we have the following stability result. 

Lemma 2.7. Let u e be a smooth solution of (PME-D) with initial data uq + e and let u be 
its uniform limit. Then u is a viscosity solution of (PME-D) with initial data uq- 



Proof. Let E be as given in (2.1) and let ip £ C 2 < l (Y,). 



1. Let us first show that u is a subsolution. First suppose that <p touches u from above 
at the point (xQ,to). We may assume that u — (p has a strict maximum at (xo,io) ^ n 
U(u) H E n {t < to} by replacing cp by 

tp(x, t) := cp(x, t) + a((x - x ) 4 — (t — t ) 2 ), a > 
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if necessary. By uniform convergence there exists a sequence (x e ,t e ) converging to (xo,io) 
such that u £ — ip has a local maximum at (x £ ,t £ ). Now if we we let 

(p{x, t) := ip{x, t) - ip(x £ ,t £ ) + u £ (x £ , t £ ) 

Then u £ — <p has a local maximum at (x e , t e ) with {u £ — (p){x £ , t £ ) = 0. We can now conclude 
by taking the limit of the viscosity subsolution property of u £ . 

2. Next we show that u is a supersolution. Let ip be a classical free-boundary subsolution 
such that tp(x, ti) ~< u(x, t±). Since the u £ 's are strictly ordered, u < u £ and hence <p(x, t±) -< 
u(x,ti) < u £ (x,t\). Now suppose (p touches u £ at some point (x2,t2), then (p(x2,t2) > 
since u £ is positive, so by continuity, there is a parabolic neighborhood of (a^^) in which 
both functions are classical and positive. By the Strong Maximum Principle, the touching 
cannot have occurred at (22,^2), a contradiction. We conclude that if < u £ so that in the 
limit ip < u £ . 

□ 



An immediate consequence of above lemma is that weak solutions are viscosity solutions 
(see Corollary 2.12). We shall introduce the precise notion of weak solutions in the next 
subsection, and summarize some results from 



2.2 Weak Solutions 



To be consistent with the setup in both [1] and [TB] , let us return to the density variable and 



consider the solution of (1.1 ) in a bounded domain with Neumann boundary condition: 

Ap m + V • GoV$) in n x R+, 

on dn x R+, 

k p(x,0) = po(x) in tt. 



(AO 



Pt 



We will see shortly that we need not worry about the fact that we are on a bounded domain, 
but for now we will let Q = f2 x R + and Qt = O x (0, t]. Following [5], 



Definition 2.8. We say p : [0, 00) — > L 1 (Q) is a weak solution of (N) if 



(i) p £ CQO^L 1 ^)) nL°°(Q t ) for all t £ (0,oo); 
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(ii) for all test functions ip G C 2,1 (Q) such that ip > in Q and dip/du = on dfl x R + , 
we have 



p(t)<p(t) = / p(0)^(0) + / / ( m + p m A^ - pV$ • Vp) 

We also define a weak subsolution (respectively supersolution) by (i) and (ii) with equality 
replaced by < (respectively >). 

/.From (J| we have existence, regularity, uniqueness and comparison principle for weak so- 
lutions: 

Theorem 2.9 (From jl]). Under the assumption that $ is C 2 in £1 

(a) The problem (N) has a unique solution; 

(b) The solution is uniformly bounded in Q and is continuous in any set £1 x [0,T]. 

(c) Suppose p{t) is a subsolution and p(t) is a supersolution, then if p Q < p i> n then 
pit) < p(t) in£lfort>0. 

The existence of solutions is obtained as the uniform limit of solutions to uniformly prob- 
lems (equicontinuity is obtained from |8j). For our purposes, a very simple approximation 
basically suffices and we summarize the relevant result in the following: 

Lemma 2.10 (From [4 ). Let u e be a solution of (N) with initial data Uq = uq + e, then 
u e is equicontinuous and there exists a subsequence which uniformly converge to u which is 
the unique weak solution to (N) with initial data uq. 

While a priori our viscosity solution is defined in all of M n , since (formally at least) solutions 
of (PME-D) should have finite speed of propagation, the boundary conditions should be 
inconsequential if we take O sufficiently large. (Later we will also establish finite propagation 



for viscosity solutions - see Corollary 2.17 ) Control on the speed of expansion of the 



support can be done via comparison with any (weak) supersolution. In particular, when $ 
is monotone (that is, when \D$>\ > except at one point where $ achieves its minimum), 
we can use the stationary profiles of the form ^(x) = (C — $)+ with sufficiently large C as 



a supersolution (see Theorem 3.1). 



Remark 2.11. Alternatively (and perhaps this is a cleaner line of reasoning), we can 
directly use the result of [11] on existence and uniqueness of solutions in all o/M n ; which 
implies in particular that the results of jll3f also apply in that setting. 



Combining Lemma 2.7 with Lemma 2.10 and the uniqueness statement in Theorem 2.9 
obtain: 



wc 
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Corollary 2.12. Any weak solution is also a viscosity solution. 



We will eventually establish uniqueness of viscosity solutions via maximum-principle type 
arguments, which culminates in the identification of the two notions of solution. 



2.3 Construction of test functions 

In this subsection we collect some test functions, i.e., (classical free boundary) sub- (super) 
solutions, to (PME-D) which will be useful for comparison purposes. In the first couple 



of lemmas (Lemmas 2.15 and 2.19), the idea is to control the $ dependence via Taylor 
expansion in a small neighborhood of a point, so that we can appropriately perturb the test 
functions for (PME) constructed in [5] and [6] for our purposes. 

The starting point is to observe that if we consider (PME-D) in some small cylinder • — 
B a (xo) x [to — a, to + a] and define vi(x,t) = a^ 1 u{a{x — xo),a(t — to)), then v% satisfies, 
in the unit cylinder -£>i(0) x [—1,1], an equation of the type 

(vi) t = (m - l)viAvx + \Dvxl 2 + b- Dv\ + 0(a)Dv 1 + a(m — l)i?iA$, 

where b = V&(xo,to). The size of the last two terms depends on the C 2 -norm of $ in Q a . 

Next we take 

v(x,t) = vi(x - bt,t) (2.4) 

Then v satisfies 

Vt = (m- l)vAv + \Dv\ 2 + 0(a)Dv + 0(a)v. (2.5) 
Next proposition illustrates the necessary perturbation one needs to perform on solutions 



of (PME) to arrive at (2.5) 



Proposition 2.13. Let u(x,t) be a viscosity subsolution of (PME) in Bi +a {0) x [—1,1]. 
Then for < a < 1, 

u\(x,t) := sup e~ at u(y,t) 

yeB a ~ at {x) 

is a subsolution of 

(PME-sub) (ui) t = (m - l)uiAui + \Dm\ 2 - a\Du\\ - au\ 

in -Bi(O) x [—1,1]. Similarly, if u(x,t) is a viscosity supersolution of (PME) in B\{0) x 
[— 1, 1], then for < a < 1, 



U2(x,t) := inf e a u(y,t) 
ytB a _ at (x) 



8 



is a su^er solution of 

(PME-super) { u 2)t = (m— l)u2^^2 + l-D^] 2 + a|-Du2| + au2 

in Bi(0) x [-1,1]. 

Proof. We only show the supersolution part. Let 112 be as given above. Suppose then that 
ip is classical and touches U2 from below at some point (xo,to). We first note that there 
exists (xi,t\) G B a _ at (xo,to) such that U2(xo,to) = e atl u(x\,ti). 

Next for any unit vector b, let us consider 

ip(x, t) = e-^-^-^ifix - (x x - x ) - ab(t -t ),t- (h - to)). 

Then we note that 1) (p{x\,t-\) = tp(xo,to) and so (u—0)(xi,ti) = and 2) by the definition 
of U2 as an infimum and by continuity of <p, in a small parabolic neighborhood of (xo,to), 
it is the case that u — (p > 112 — </5>0;we therefore conclude that (p touches u from below 
at and so we have, 

[^A^ + |V^| 2 ](x ,t ) = e Q *i[^A^+|V^| 2 ](a;i,ti) 

< e at ^ t {xiM) 

= [ip t ~ aip - ab ■ Vip](xo,to). 

Now the desired inequality is achieved by setting b = j^r(^O) *o)- 

Indeed the above calculation shows that if u is a viscosity supersolution of (PME), then U2 
should be a viscosity supersolution of (PME-super): If a classical free boundary subsolution 
(p of (PME-super) crosses U2 from below, then the corresponding (p is a subsolution of 
(PME) and crosses u, yielding a contradiction (there is no distinction between the interior 
and boundary cases). □ 

Lemma 2.14. [CVJ Consider the function 

H(x,t;A,u) = A(\x\ +ujt- B)+ 

with R/2 < B < R. Then u is a classical free boundary supersolution of (PME) in the 
domain {\x\ < R} x [^(R - B), 0] if 

uj , w , R — B 

_ > l + 2 ( m -l)(d-l)_ s _. 

Proposition |2.13 and Lemma 2.14| yields the following: 
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Corollary 2.15. Let us fix xq G W 1 and let H be given as in Lemma 2.1^. Then the inf 
convolution of H, given as 

H(x,t;a) =e at inf H(y,t) 
y&B a - at {x) 



is a classical (free boundary) supersolution of ( |2.5| ). Consequently, there exists C = Co 
which only depends on the C 2 -norm of in B±(xo) such that 

H(x, t) := aLB[oT l {x — x$ + b(t — to)), a~ (t — to); Co) 

is a classical (free boundary) supersolution of (PME-D) in Q a := B a {xo) x [to — a, to]. 



Proof. By Lemma 2.6 and Proposition 2.13, H_ is a viscosity supersolution of (PME-super), 
so it is sufficient to show that it has the required regularity. For this, note that for simplicity 
we have only taken the supremum over space and the reader can readily check that in this 
case the infimum for H(x, t; a) is achieved at the point y which minimizes |y| subject to the 
constraint that \y — x\ = a — at, and thus an explicit expression is possible for LI. □ 

Remark 2.16. In fact due to the explicit form of H it follows that the free boundary velocity 
of H is given by 

-> DH 
V = uj + b- —— + Ca. 
\L)ti 



By comparison with these supersolutions, we immediately obtain 

Corollary 2.17. Any viscosity solution has finite propagation speed and is bounded in a 
big ball in any local time interval. Further, if $ is convex, then via comparison with the 



stationary solutions of the form given in Theorem 3.1, the above holds globally in time. 



For the analysis in section 3, we will make use of the Barenblatt profiles (see e.g., [CV] and 
[V])- 

Lemma 2.18. [Barenblatt] Let B(x,t;r,C) be the family of functions 

(C(t + T) 2X - K\x\ 2 )+ 



B(x,t;r,C) ■■ 
with constants A, K, C, r > such that 



t + r 



A = ((m - l)d + 2) _1 , 2K = X. 
Then B(x,t;r,C) is a classical (free boundary) solution of (PME). 



Using Proposition 2.13 (see also the proof of Corollary 2.15) once again, we obtain the 
following: 



10 



Lemma 2.19. Let us fix x$ £ W 1 and let B be a Barenblatt function. Then there exists C 
which only depends on the C 2 -norm of & in Bi(xq) such that 

ip{x, t) = ae~ Ca ^ ta) sup B(a~ l {y-XQ + ab{t-t )),a~ l {t-tQ)) 

y£B Ca _ c(t _ to) (x) 

is a classical (free boundary) subsolution of (PME-D) in Q a := B a {xo) x [to — ct,£o]- 
Remark 2.20. The reason for taking the hyperbolic scaling is because we will have occasion 



to require rather fine control on the boundary velocity (see Lemma 2.24) an d this is the 
scaling which preserves the velocity - in contrast to the parabolic scaling, which dramatically 
reduces the effect of the drift <]? in the bulk (the positivity set), but unfortunately at the cost 
of severely disrupting the boundary velocity. 

To establish the Comparison Principle, we will need the following weak analogue of (ii) 
in the definition of viscosity supersolutions for subsolutions, the proof of which utilizes an 
approximation lemma from [5] . 

Lemma 2.21. Let u be a viscosity subsolution of {PME-D), and let ip be a classical free 



boundary supersolution from Lemma 2.15 which lies above u at some time to- Then (p 



cannot cross u from above at a later time t > t$. 

Proof. Let ip and u be as described in the statement, and suppose that ip touches u from 



above at some point (xq, to). From Lemma 2.15 we have that <p is given as the inf convolu- 



tion of some spherical traveling waves from Lemma |2.14| which we denote ip. Further, let 
us suppose the infimum is achieve at (xi,t\) so that 1) (p(xo,to) = ip(xi,ti) and 2) by the 
definition of <p as an inf convolution, the translated function 

■tp(x, t) = ^{x + (xi - x ), t + (ii - t )) 

also touches u from above at the point (xo,^o)- From Lemma 4.4 in [5], we know that 
ip can be given as the monotone limit of classical positive supersolutions, and hence the 
same is true of ip: I.e., there exsits ip e \ ip, with ip e > classical. But since u cannot 
touch tp e by the Strong Maximum Principle, we obtain in the limit that u < tp, which is a 
contradiction. □ 



2.4 Comparison Principle and Identification with Weak Solution 

Here the outline of the proof closely follows that of the corresponding result for (PME) 
(Theorem 2.1 in j^]): We will give an abridged version of the proof, pointing out main steps 
and modifications for our problem. 

Theorem 2.22. [Comparison Principle] If u is a viscosity subsolution and v is a viscosity 



supersolution in the sense of Definition 2.4 with strictly separated initial data, uq -< vq, then 
u(x,t) < v(x,t) for every (x,t) £ Q. 
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Proof. 1) [Sup and Inf Functions] For given 5 > and r > small with r > 5, we introduce 
the regularized functions 

W(x,t) = _ inf v(y,r) 
B r _ st {x,t) 

and 

Z(x,t) = sup u(u,t) 

B r (x,t) 

First note that and Z preserve properties of v and it: 

o W is a supersolution and Z is a subsolution; 
o Z(-,r) -< W(-,r) for r sufficiently small. 



For a proof of the first item see [5], Lemma 7.1 or the proof of Lemma 2.15). The proof of 
the second item relies on the fact that 



o The support of a viscosity subsolutions and supers olutions evolve in a continuous way. 
Here continuity is understood as continuity in the Hausdorff distance (in time) of the 
positivity set. 

The proof of this can be done by comparison with the supersolutions (respectively subso- 



lutions) constructed in Lemma 2.15 (respectively Lemma 2.19). We omit the details since 



with replacement of barriers it is no different from the proof of Proposition 6.2 in [5]. 

Thus if we can prove that W stays above Z for all choices of r and ^(sufficiently small), then 
we may take 5 — > and then r — > to recover the conclusion for it and v. First let us note 
that due to the Strong Maximum Principle, W cannot touch Z from above, and therefore 
we are reduced to the analysis of a first contact point of W and Z at some Pq = (A, to). 

The key usefulness of Z and W lies in the fact that they enjoy interior /exterior ball prop- 
erties: 



o The positivity set of Z has the interior ball property with radius r at every point of its 
boundary and at the points of the boundary of the support of u where these balls are 
centered we have an exterior ball; 

o The positivity set ofW has the exterior ball property with radius less than r — 5t (since 
in this case we really have an exterior ellipsoid in space-time) and at the points of the 
boundary of the support of v where these balls are centered we have an interior ball. 

For detailed proofs of these statements we again refer the reader to [5] . 
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Figure 1: The geometry at the contact point 



2) [The Contact Point] The first contact point Pq = {xq, to) is located at the free boundary of 
both functions. Therefore by the definitions of Z and W, there is a point Pi = (xi,ti) on the 
free boundary of u located at distance r from Pq and there is another point P2 = (£2^2) 
on the free boundary of v at distance tq = r — Sto from Pq. Let us also denote by Hz 
(respectively H\y) the tangent hyperplane to the free boundary of Z (respectively W) at 
Pq. (see Figure [T]) 

Lemma 2.23. Neither Hz nor Hyy is horizontal. In particular, one can denote the space- 
time normal vector to Hz, in the direction of P\ — Pq, as (e n ,m) G W 1 x R where \e n \ = 1 
and —00 < m < 00. 

Proof. It is enough to show that t\ > tQ — r (i.e., T(Z) cannot propagate with infinite speed) 
and t% < to + r (i.e., T(W) cannot propagate with negative infinite speed). The desired 
conclusion then follows by the ordering of Z and W. 

We first show that ti < to + r - Otherwise H\y is horizontal and after translation we have 
Pq = (0, — r) and P2 = (0,0). Moreover Q(v) has an interior ball at P2 with horizontal 
tangency with radius < r' < r. Now in any parabolic cylinder 

C v = {(x, t) : \x\ < 77, -r] 2 < t < 0} 
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A t 




i = t 2 



P(V) 



lb • 



v > o 



Figure 2: Infinite, negative speed 



with bottom edge contained in the interior ball (which can be achieved by taking -n < r'), 
we have by continuity that v > M > on that edge, (see Figure [2| 

On the lateral boundary of it may be the case that v = 0, so we will have to compare 
with a subsolution with support strictly contained in —rj < \x\ < r\ at time t = —rj 2, and 
still contains in its support at time t = 0, which rigorously implies that v cannot contract 
sufficiently fast for (0, 0) to be a free boundary point. The necessary subsolution can be 



constructed as the one in Lemma 2.19 adjusted for the parabolic scaling. 



The case t± > to — r follows mutatis mutantis from the arguments in [6] and [5] (see [6], 



Lemma 4.2 or [5], Lemma 8.2) with the barrier in Lemma 2.15 replacing the barriers used 
therein. □ 



3) [Non-tangential Estimate] The next lemma states that the normal velocity V of T(Z) at 
(xo,to) satisfies, in the viscosity sense, 

DZ 

V [xoM) < (\DZ\+D<t>(x ) • 7^r)(*o,to). 
Lemma 2.24. Let x n := x ■ e n , and consider a non-tangential cone K := {x : x n > 
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k\x\ with k > 0}. Then we have 

Z(x + x,t ) , , 

lim > m — D<P(xo) ■ e n . 

x£K, x^O X n 

Proof. The argument is parallel to the proof of Lemma 4.3. in [6]; the only difference for 
us is taking into account the change of reference frame introduced by the drift given by <£. 



This is ensured by the local nature of the construction of our barrier in Lemma 2.15| which 
replaces the corresponding barriers used in [6]. 

□ 



4) [Conclusion] Due to Lemma 2.24, we may place a small subsolution <p from Lemma 2.19 



below Z at Pq with speed close to m (again see Remark 2.20, which assures us that our 
subsolutions are constructed so that this is possible) such that it crosses anything with 
speed m! < m. Since </? is also below W and hence v (after a small translation), v must 
expand by at least m', but then r(W A ) has speed m! + 5 > m at Po, yielding a contradiction 
to the fact that Z touched W from below at Pq- D 



We can now establish uniqueness of viscosity solutions: 

Theorem 2.25. The problem (PME-D) admits a unique solution in the class of viscosity 
solutions as defined in Definition \2.J\ for continuous and nonnegative initial data. This 
solution coincides with the continuous weak solution. 



Proof. The existence of a continuous weak solution can be provided as the uniform limit of 
classical solutions with initial data uq E = uq + e, and by Lemma 2.7 such a limit, which 



we will denote by U, is also a continuous viscosity solution. Further, by comparison with 
uo,e and taking a limit, it is clear that such a limit U is also a maximal viscosity solution. 

Uniqueness would follow if we can show that any other viscosity solution u also cannot be 
smaller than U. For this purpose, consider u n (x,t) with initial data u n (x,0) := (no — -)+• 
Now consider positive u^ 1 such that \u £ ^ — u n \ < - in W 1 x [0,T]. It follows from Lemma 



2.10 that u^ 1 uniformly converges to U2(x,t), which is then a continuous weak solution of 



(PME-D). Therefore, by uniqueness of weak solutions, U2 is equal to U. On the other hand 
by Theorem 2.22 u n -< u and thus U = XJ% < u. Hence we conclude. □ 



Using Theorem 2.22| and Theorem 2.25 we can in fact prove a stronger comparison theorem 
for viscosity solutions (see [5], Theorem 10.2): 

Theorem 2.26. Let u\ and U2 be respectively a viscosity subsolution and a viscosity super- 
solution of (PME-D) in some parabolic cylinder Q with initial data uq,i and uo,2 such that 
Uo,i(x) < itQ^- Then u±(x,t) < U2(x,t) for all (x,t) G Q. 
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This comparison theorem in particular allows us to restrict attention to only the (classical 



free boundary) supersolutions used to establish Theorem 2.22 , and consequently, we can now 



strength Lemma 2.21 to enable comparison with any classical free boundary supersolution: 



Lemma 2.27. Let u be a viscosity subsolution of (PME-D), and let tp be a classical free 
boundary supersolution which lies above u at some time to. Then <p cannot cross u from 
above at a later time t > to. 

Proof. Let us first replace <p by a viscosity solution of (PME-D) with the same initial data 



which we denote v. Since tp is a viscosity supersolution by Lemma 2.6 we have that p > v 



by Theorem 2.26 Finally, u < v by Theorem 2.22 □ 



3 Convergence to Equilibrium 

We begin by discussing the set of equilibrium solutions to (PME-D) and reviewing some 



known results. Since by Theorem 2.25 the unique viscosity solution coincides with the 



continuous weak solution, we may carry out our discussion in the context of weak solutions. 

The set of equilibrium solutions and uniform convergence of solutions to the equilibrium 
are established in [4\. Below we state the corresponding result in the pressure variable. 

Theorem 3.1. [Theorem 5.1, US] The set of equilibrium solutions for (PME-D) is given 
by 

S = {W e C{R n ) : W > in O, and for every x£(l, either W{x) = 
or W + & = C for some constant C in a neighborhood of x} 

Further, given uq, there exists a unique W{x) G S such that u(x,t) uniformly converges to 
W as t — > oo. 

It is fairly immediate that S is contained in the set of equilibrium solutions; the converse 
containment and the convergence statement are established based on a L 1 contraction result 
in @]. 

Under the assumption that is convex. Note that the density function p{-,t) given in 



(1.2) preserves its L 1 norm over time. Therefore there is a unique equilibrium solution 



Uoo = (Co — 3>)+ to which u(-,t) uniformly converges as t — > oo, i.e., the one with 

(u ) 1/{m - 1 \x)dx= [ ( Uoo )^ m -\x)dx. 



An explicit exponential rate of convergence is then derived in [13] by the entropy method: 
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Theorem 3.2. Suppose that & is strictly convex, i.e., there exists a constant ko > such 
that x ■ [(Hess$>(x))x] > ko\x\ 2 for x 6 W 1 . Let Uqo = (Co — <£)+ be the equilibrium 
solution to which our solution u(x,t) converges as t — > oo. Then there exist constants 
K, a > depending on m, ko and the L 1 norm of uq such that 



\u(x,t) - Uoo{x)\dx < Ke at . 
Remark 3.3. In fact the estimate in I13f is given in terms of the pressure variable p. 



Due to Corollary 2.17, for a convex (and in fact monotone) potential <E> u is uniformly 
bounded with its support contained in a compact set for all times. This allows us to derive 
the estimate for u from that of p for 1 < m < oo. Further, due to the equivalence of all 
LP norms in our setting, we will take some liberties in passing between u and p in our 
estimates. 



3.1 Convex potential 

As an application of the viscosity solutions theory, we will convert the L 1 estimate in 



Theorem 3.2 into a pointwise estimate (see Lemmas 3.4 and 3.5); such an estimate in turn 



will yield a quantitative estimate on the rate of the free boundary convergence (see Theorem 



3.6). 



Rescaling u by cu(x, ct) if necessary, let us assume for the rest of this subsection that u < 1 



and maxA$ < 1 on our domain of consideration, which is bounded (see Corollary 2.17) 
and we assume to be {|x| < R} for some R > 0. 

Lemma 3.4. [Uniformly Bounded From Below] There exists a sufficiently small constant 
k > 0, depending only on m and n, such that the following is true: 

Suppose, for (xo, to) G M n x (0, oo) and for < a < 1, 

u(-,to)dx > a k . 
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Then u(-,to + a) > a k ' in B a (xo). 
Proof. 1. Let us define 

u(x, t) := u(a(x — xq), a 2 (t — to))- 

Then u is a supersolution of 

^ = (m - 1)V>At/! + \DiP\ 2 - Ca(\D^j\ + V), (3.1) 
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where C is a constant depending on the C 2 -norm of (near xq). Below we will construct 



a subsolution of (3.1) to compare with u in order to establish the lemma. 



2. Let us consider w(x,t) which satisfies 

w t = (m-l)wAw + \Dw\ 2 -Ca (3.2) 

in the weak sense (see e.g., [9]), with initial condition w(x, 0) := (1 — a)u(x, 0)x\x\<ii where 
rh — 1 = (1 + a){m — 1). 

Then, say for t > 1/2, w(-,t) is Holder continuous due to [9]. Since u is bounded by 1, 
so is w(x,0), by the Maximum Principle. Further, note that any solution of the (PME) is 



now automatically a supersolution of (3.2) and therefore, using an appropriate Barenblatt 



profile as a supersolution of (3.2), one can check that 

n 1/2 (w) C {|x| < 2}. (3.3) 

Moreover, integration by parts yields that 

w(x,t)dx = / w(x, 0)dx — at, 



and in particular we deduce that e.g., f w( x, l/ 2)dx > a k /2. Let x* be the point where 
ii;(-,l/2) assumes its maximum, then from (3.3) we see that w(x*,l/2) > C n a k for some 



dimensional constant C n . Due to the Holder regularity of w, 

w(-,l/2)>(C n /2)a k in B ak2 (x*). 
where k 2 = j~ 1 k where < 7 < 1 depends only on m and n. 

Let U{x,t) := B{x - x*,t;T,C), where B{x,t;r,C) = Mt+r)^-K\x\*) + ig thg Barenblatt 
profile given in Lemma |2.18[ with rh as the permeability constant and the conditions 

Ct 2X - x = a k /A (height), and ^JC/Kt x < a k2 /2 (the size of initial support). 

Note that, in its positive set, U is concave. Therefore U(x,t) := (U(x,t) — Cat) + satisfies, 
in its positive set, 

U t = Ut-Ca 

= (m-l)(U- Cat)AU + (m - l)atAU + \DU\ 2 - Ca 
< (m-l)UAU + \DU\ 2 - Ca. 



Therefore U(x,t) is a subsolution of (3.2) for positive t. Comparison with U(x,t) and 
w(x,t + |) yields that 

w^a' 1 ) > a k ' in ^(0), 
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if k' is chosen sufficiently large. 



3. Observe that v := (1 + a)w satisfies 

v t = < (m - l)vAv + j^\Dv\ 2 - Ca(l + a) 

< (m - l)uAw + \Dv\ 2 -Ca{\Dv\ 2 + 1) 

< (m-l)vAw + |£>u| 2 -Ca(|Dv| + l) 



Hence u is a subsolution of (3.1), and further v(x,0) = (1 — a 2 )u(x,0) < u(x, 0). Therefore 
by the comparison principle (for weak solutions, see |9j) we have v < u. From previous 
discussions we obtain that u(-,o _1 ) > a k in -Bi(O), and thus 



u(-,a) > a k> in B a (0), 



as desired. 



□ 



Next lemma establishes a uniform upper bound on u in terms of its L 1 norm. Here we would 



address our equation in pressure form, i.e., (1.1), to invoke regularity theory for divergence 
form operators studied in |15| . 

Lemma 3.5. Let K be a compact subset ofW 1 with u = outside of K for all time. Then 
there exists a constant C > depending on m > 1, sup p and maxA$ := M\ such that the 
following holds: 

If 

p(-,t)dx < c for ti < t < t 2 := t\ + log(l/c ). 



Ib, 



'B o (0) 

then p{;t 2 )< Cc Q /n+1 in B^O). 

Proof. 1. We proceed by induction. Suppose that u < a = c^ n+1 2 k in -Bc(O) x [O,^] where 
k > is chosen such that 

4Ccl /n+1 < a < 1, (3.4) 
with C to be determined later. Our goal is to show that 

u < a/2 in B c(1 _ a/2) (0) x [l,t 2 ]. 

Then the desired result is obtained by iteration, beginning with a = 1 and continuing until 



a reaches the lower bound in (3.4). Note that the total number of iteration for this process, 



therefore the total time we need for the desired result, is of order log(l/co). 



19 



2. Let pi(x,t) solve our equation (1.1) with initial data po + a/10 in S := Bc(0) x [0, oo) 



with boundary data corresponding to it + a/10. Observe that f(t) = yAae* is a supersolution 



of (1.1), since maxA$ = 1 due to our normalization. Thus by the comparison principle 
pi ranges from a/10 to 4a for < t < 1. Therefore, treating the diffusion coefficients as a 
priori given, p(x,t) := a~ 1 pi(ax, t) solves a quasi-linear equation of divergence form with 
diffusion coefficient of unit size: 

(P) p t = V • (b(x,t)Vp + pV$), where b(x,t) = p e [1/10,4]. 

In particular, we can decompose p := pi + p2 where pi solves (P) (which is linear) with 
initial data po/a + 1/10 and boundary data zero, and p 2 solves (P) with initial data zero and 
boundary data 4. We claim that both of them stays smaller than 1/4 in -8(7/0-1/4(0) x [0> !]• 

For pi we have 

/ pi(-,t)dx< pi(;0)dx = c a~ (n+1) for all t > 0, 



Due to (3.4) and the Holder regularity of p~\ (see [15]) we have 

pi(; 1) < Cc a-(™ +1 ) < 1/4 in B 1/a _!(0). 



if we choose C > sufficiently large in (3.4) corresponding to the Holder regularity for 
solutions of (P) for t > 1. 

As for p~2, arguments with test functions in the weak formulation of (P) (not very different 
from the case of the heat equation) in combination with the Holder regularity estimates 
yield that, for sufficiently large C > 0, 

p 2 (x,t) <4te" 3 /* in B c/a _ c/ M 

with a dimensional constant C. In particular p 2 {x, 1) < 1/4 in -Bc/a-C/4(0)- 

Hence we obtain p(x, 1) < 1/2 in B c / a _ c / 4 (0). Scaling back to our original density function 
p(x,t), we conclude that 

p(.,l) <o/2 m£ c(1 _ a/4) (0). 

Due to our assumption on po, one can go through the above argument starting at any time 
t = t £ [0,^] instead of t = 0, obtaining that p(-,t) < a/2 in Bcn-a /2) (0) x [1, ^2] • 

5. Repeating the above argument with a 2 '■= a/2 starting at t = 2, we get 

p < a 2 /2 = a/4 in B C (i-a/4-a/8)(0) x [3,*2]. 



If we iterate up to of order logco times, then a reaches the lower bound in (3.4), and we 
arrive at the desired result. 

□ 
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Next we use the uniform bounds obtained above to investigate the rate of free boundary 
convergence. 



Theorem 3.6. Let $ and Uoo be as in Theorem 3.2. Then forT > sufficiently large, T(u) 
is in the Ke~ a2t -neighborhood o/r(u t>0 ). Here 1/K,a 2 > is a sufficiently small constant 
depending on m, sup«0; ko, Mi, A := m.in^^ >Co \D$\ and n. 

Proof. 0. We first show that Tt(u) is close to the equilibrium profile from the inside, i.e., 

{x : d(x, R n - n(woo)) > Ce~ ct } C n t (u). (3.5) 

After t = T > Clna, L x -average of u in B a (xo) is bigger than Ca k when the center xq is in 
S^a := < C - a k }. Hence Lemma |3~4| yields that u(x ,T) > C 2 a k ' with T = a k / 2 . As a 



result we conclude that T(u) lies outside of tt a after t = Clna. 

It remains to show that is also close to the equilibrium profile from the outside. This 

will be more involved. 



1. Let D a := {$ > C + a} and let b(t) := e~ at . Due to Theorem [3^2| we have 

p{x,t)dx < Ke~ aT for t > T. 
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Take any point xq such that B Cb ^(xo) C D^t) if T is chosen large. Hence one can apply 
Lemma 13.51 to 

t) := u(b(T)(x - x ),b(T) 2 t + T) 
with t\ = and cq = Ke~ aT to obtain 

n(-,aT + /3) < K' e - aT/{n+1) in Bi(0). 

Since xo was arbitrarily chosen, repeatedly using Lemma |3.5| and scaling back we get 

u(x,t) < Ke~ aT/in+1) in D 2b{T) x [Ti,oo),T := (1 + Ca)T. (3.6) 



2. We claim that estimate (3.6) yields that 

T t (u) C M n - D a[t) with a{t) := Ke"" 2 ' for t > T 



(3.7) 



2C 



where ot2 = l/i^ 3 and K > — — , where Co is given in Corollary 



2.15 



We prove the claim 

by induction. Due to Corollary 2.17 the claim is true up to a sufficiently large time t = Tq 
if K is sufficiently large (in particular a 2 < a). Next let 



T* := sup{t : (3.7) holds for < t < t } > Tq. 
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We want to show that T* = oo. Let us choose xq G dD a , where o = Ke~ a ' 2T * . Since 
is convex with \D<&\ > in D a , there is an exterior ball B\(xi) outside of D a such that 
xo G dBi(xi). Due to (3.6) and due to the fact T* > To, we have for a := 1/K, 

u < a/K in £ := [B 1+ - a { Xl ) n B^xq)] x [T*, oo). 

Let us first make a heuristic argument. Suppose that xq G Tt*(u) (by the definition of T* 
such an x$ exists). Then at xq we have DQ(xo,T*) pointing in the direction of x\ — xo, 
which is the inward normal of T(u) at (xq,T*), which is parallel to — jjj^(xo,T*). Hence 
formally the normal velocity of T(u) at (xo,T*) should be 

V = (\Du\ - \D<S>\)(x ,T*) < 0(a) - A < -A/2 (3.8) 

if K > 2/ A. This way one can conclude that the furthest point of Tt(u) from Dq shrinks 
away with the normal velocity less than — A/2 at t = T*. This would yield a contradiction 
to the definition of T* since 

-(Ke~ a2t Y = Ka 2 a < A/2. 

(see Figure [3]) 

The rest of the proof consists of a barrier argument to establish an appropriate version of 



3. Consider H from Lemma with u = 8/K, A = A/K, B = 3/4, R = 1 and define 
with a = d and (xo,^o) = {xo,T*). Then H is bigger than u on the 



2.15 



H as in Lemma 

parabolic boundary of E n {T* < t <T* + e} for sufficiently small e > 0. We then compare 
H with u in £ n [T* , T* + e\. Observe that the spatial outward normal of T(H) at (xq, T*), 
which is —DH/\DH\, points toward x\ — xq. Hence the outward normal velocity of T(H) 
satisfies (see Remark 2.16) 



V = uj- \D$\(x ,T*) + Cd< -A/2. 

if K > 2C/A. Since u < H in S n [T*, T* + e), we conclude that 

u(-,T* +e) = in B Ae/2 (x ). 

Since xq G dD a is chosen arbitrarily, we conclude that u(-,T* + e) vanishes in Ae/2- 
neighborhood of D a , which includes D a ^ +£ ^ if K > 2/ A for any e > 0. 



This contradicts the definition of T*. 



□ 
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Figure 3: The construction of barriers 
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3.2 Monotone potential 

Suppose $ is monotone, i.e., 

\D$>\ > except at x = xq where $ obtains its minimum. 

We are not able to yield quantitative estimates on the rate of free boundary convergence, 
due to the lack of available L 1 -estimates. However we state the theorem below to illustrate 
that if we neglect the rate of convergence, then considerably simpler arguments already 
yield the free boundary convergence of u as t — > oo, 

Theorem 3.7. Suppose $(ac) satisfies \D$\(x) > except at x = 0, where <3? achieves 
its minimum. Let Uoo(x) := (Co — <&)+(a?) where Co is chosen such that ||woo|Ili/™-i = 
||wo|lx,i/m-i. Then T t (u) uniformly converges to r(« 00 ) in the Hausdorff distance, as t — > oo. 



Proof. 0. We first verify that for any compact subset K of O(tioo) there exists T > such 
that Tt(u) lies outside of K for t > T. This is immediate from the uniform convergence of 

U to Uqo. 

1. It remains to establish the statement from the outside of r2(woo). 

Since $ is monotone the sets {x : $(x) < r} starts from a point (when r = min<I>) and 
extends all the way to R n . In particular, there is only one equilibrium solution to which 
u converges as t — > oo. Choose T > 0. Recall that, by the comparison principle, the set 
£l(u) n [0, T] is compact in W 1 x [0, oo). Therefore there exists a largest C(T) > such that 
{x : <J?(x) < C} contains Ot("u)- We will show that fJ(it), after a while, eventually shrinks 
away so that after T\ > T we have 

Oti(«) C {x : $(s,t) < C 2 } with C 2 < C. (3.9) 

Let us choose x$ G d{Q(x) < C}. Due to the regularity of $ there is a ball £? r (:ci) C 
{<&(x) > C} such that xo G dB r (xx), where r > is independent of the choice of a^o- Due 
to [3], u(x,t) uniformly converges to zero in {x : &(x) > C} with C > Cq. Therefore there 
exists To such that 

u < r/2min |-D$|(x) in S x [T , oo), where E := {x : C < C 2 < < C + 1}. 
(It is noted that the condition is to ensure that, formally, \Du\ < \D$>\ at (xq,Tq)). We now 



argue as in step 3 of the proof of Theorem 3.6 to construct a radially symmetric barrier 



supported in (B2 r (xi) — B r n\(xi)) X [To, To + 1] with r'(t) < 0, boundary data zero on 
dB r n\(xi) and 2r on dB2 r (xi), to demonstrate that at T(u) stays out of B r / 2 (xo) after 



t = Tq + 1. Since a?i is an arbitrary point of d{&(x) < C}, we arrive at (3.9). 



□ 
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